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Polarization of asteroid (387) Aquitania: the newest member of a 
class of large inversion angle asteroids 

Joseph Masiero^, Alberto Cellino^ 
ABSTRACT 

We present new imaging polarimetric observations of two Main Belt aster- 
oids, (234) Barbara and (387) Aquitania, taken in the first half of 2008 using the 
Dual-Beam Imaging Polarimeter on the University of Hawaii 2.2 meter telescope, 
located on Mauna Kea, Hawaii. Barbara had been previously shown to exhibit a 



very unusual polarization-phase curve by lCellino et al.l (120061 ). Our observations 
confirm this result and add Aquitania to the growing class of large inversion 
angle objects. Interestingly, these asteroids show spinel features in their IR spec- 
tra suggesting a mineralogical origin to the phase angle- dependent polarimetric 
features. As spinel is associated with calcium-aluminum-rich inclusions and car- 
bonaceous chondrites, these large inversion angle asteroids may represent some 
of the oldest surfaces in the solar system. Circular as well as linear polarization 
measurements were obtained but circular polarization was not detected. 



1. Introduction 

The radiation we receive from asteroids at visible wavelengths is sunlight scattered by 
the solid surfaces of the objects. The scattering process polarizes the emerging photon 
flux, with the most general state of polarization being partial elliptical polarization. In 
the case of asteroid scattering, linear polarization dominates over circular and is modulated 
by the properties of the surface (e.g. albedo, texture, composition, regolith size) and the 
illumination conditions. Measuring the degree of linear polarization can diagnose physical 
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conditions of the scattering surface and is complementary to photometry and spectroscopy 
for the remote analysis of small solar system objects. 

A common way to quantify the polarization of light is by using the Stokes vectors J, 
Q, U, & V, where I is the total intensity, Q the amount of light polarized in the 0° or 
90° planes, U the light in the 45° or —45° planes, and V the light polarized circularly with 
left- or right-handedness. In particular. Stokes Q and U fully describe the state of linear 
polarization. For cases of linear polarization, Q and U can be converted into a measure of the 

total degree of polarization (P) and the position angle of polarization (9) using: P = — — j 

and 9 = | arctantZ/Q. 

The reference for the position angle of polarization is typically celestial north, though 
definitions in the instrument frame or the galactic frame can also be used if convenient. In 
asteroid polarimetry it is found that with few exceptions the orientation of the plane of linear 
polarization is perpendicular or parallel to the scattering plane: the plane encompassing 
the Sun, the target, and the observer. Since most scattering cases result in light polarized 
perpendicular to the scattering plane we define that direction as the reference for polarization 
position angle. The new degree of polarization is then defined as 

where I± and /y are the intensities of light perpendicular and parallel to the scattering 
plane, respectively. For this definition, Pj. will be positive for the "normal" scattering case, 
or negative if the scattered light is polarized in the plane of scattering. This allows us to 
refer to polarization of an asteroid as either "positive" or "negative" and still be completely 
descriptive. 

The primary observable in asteroid polarimetric studies is the change in percent polar- 
ization Pr as a function of phase angle (the angle between the Sun and the Earth as seen from 
the asteroid). Close to zero phase angle Pr tends to zero, but as the phase angle increases Pr 
becomes increasingly negative. This is the so-called branch of negative polarization. After 
reaching a maximum negative value usually between 8° and 10° phase angle, Pr decreases 
in absolute value (that is, becomes depolarized) and reaches zero at the so-called inversion 
angle (oq) usually between 15° and 20° phase. The maximum value of negative polarization 
is an important parameter, usually indicated as Pmin, and varie s mostly between —0.5% and 



—2.0%. Measurements of Pmin can be used to calculate albedo (IZellner &: Gradidll976l). and 



can b e used along with inversion angle as an indication of surface texture (jPoUfus fc Zellner 



1979|). 



Beyond the inversion angle Pr becomes re-polarized in the positive sense, increasing 
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nearly linearly with phase. The slope [h) of the branch of positive polarization for phase 
angles greater than the inversion angle is another impo rtant parameter, si nce it is known 



to be diagnostic of the albedo of the surface (see, e.g. ICellino et al.l Il999l . and references 
therein). The maximum value of positive polarization is reached at values of phase which 
are well beyond the maximum values attainable by main belt asteroids observed from Earth. 
Only near-Earth objects are occasionally visible over very large intervals of phase angle, 
and the maximum value of positive polarization is observed in some cases to occur at phase 
angles of the order of 80 - 100 degrees. 



Shkuratov et al.l (119941 ) review the many possible physical models that have been used to 
try to explain the negative polarization seen in asteroid phase curv e s. Th e most promising are 
the coherent backscattering models, described in e.g. iMuinonenI (119891 ). which can explain 
both the photometric opposition surge observed for atmosphereless bodies as well as the 
negative polarization seen for asteroids at small phase angles. The polarimetric effect in this 
theory has a strong dependence on the size and spacing of the scattering particles meaning 
that surface chemistry and mineralog y should play a la r ge rol e in determining the polarization 
as a function of phase. According to iMuinonen et al.l ((20021), coherent backscattering is the 
key physical process generating both these effects. 

As part of an extensiv e campaign to characterize the polarization-phase curves of over 
100 asteroids ICellino et al.l (120051 . l2006l ) found an interesting case in the asteroid (234) Bar- 
bara, which was found to exhibit an unusual polarimetric behavior. Whereas most asteroids 
have an inversion angle of around 20° Barbara shows strong negative polarization (~ — 1%) 
for phase angles larger than 25°. The authors proposed that this strange behavior was due 
to unusua l surface properties, hke ly the surface mineralogy driving the rare Ld-type clas- 
sification. iGil-Hutton et al.l (120081 ) performed a follow-up investigation of other L/Ld-type 
asteroids, finding four that show these unusual pol arization properties. While not all L-type 
asteroids show this ef fect (e.g. (12) Vi c toria, see ICellino et al.l (l2005l )). both Barbara and 
(980) Anac qstia from iGil-Hutton et al.l (120081) 's survey show strong spinel features in their 



IR spectra (ISunshine et al.ll2007l ). iBurbine et al.l (119921 ) identified strong spinel features in 
the IR spectra of both (387) Aquitania and Anacostia, so we investigated Aquitania to look 
for the same strange polarization signatures seen in other spinel-rich L-types. 



2. Observations 

Observations of the asteroids (387) Aquitania and (234) Barbara were conducted on 
four nights spread from January to June 2008 on the University of Hawaii 2.2 m telescope 
located on Mauna Kea, Hawaii. Polarizations were measured with the Dual-Beam Imaging 
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Polarimeter (DBIP. iMasiero et al.l (120071 )). a broad-band CCD imaging system sensitive to 
both linear and circular polarizations simultaneously. DBIP has very low instrumental sys- 
tematics and can reach polarization precisions of better than 0.1%. DBIP uses an IR-blocked 
clear filter which transmits from 400 — 700 nm, a close approximation to a Sloan q' + r' filter . 
Although there is some evidence for color dependence of polarization (jCellino et al.ll2005l ). 
the difference in measured polarization between Vy and is usually within error. V and 
R band polarizations of Barbara showing this consistency can be seen in Fig [H Using a 
broader filter does not impair our measurements of the bulk polarization properties. 

Stokes parameters were measured for our targets using a beam-swapping pattern of 
waveplate alignments. Each image measured complimentary values of a single Stokes vector 
(i.e. I + Q in the north beam and I — Q m the south beam; I + U north and I — U south; or 
I + V north and I — V south) which were then swapped in a second image. By combining 
both images we were able to completely remove time-dependent and flat-field effects on the 
measured polarization values. Individual exposure times were adjusted so that each image 
of our target had a peak count value on the CCD between 20, 000 and 40, 000 counts per 
pixel to provide enough photons for good statistics while avoiding non- l inear issues near 
saturation. Image reduction was completed as described in lMasiero et al.l (120071 ). 



DBIP's native polarization measurement is of the fractional values Q/I = q and U/I = 
u, which can be converted to P and 6, as well a.sV/I = v. Table [1] presents the asteroid 
name, UT date of observation, apparent V magnitude, total exposure time, phase angle (a), 
measured linear polarization [Pr), angle of polarization (6) and measured circular polariza- 
tion. Total hsted exposure time includes all six waveplate positions needed to develop a full 
measurement of q, u, and v. The linear polarization measurements and angles of polarization 
have been rotated into the frame of the asteroid scattering plane. The circular polarization 
values in all cases are consistent with a zero signal to within 2.5a, and thus we do not detect 
any circular polarization of the light scattered from these asteroids. Circular polarization 
is predicted to appear on asteroids that have surfaces containing powd ered metals, espe- 
cially those with highly irregular shapes (iDegtjarev fc Kolokoloval Il992[). however circul ar 



polarization as of yet has not been detected from an asteroid (e.g., iMuinonen et al.ll2002l ). 



Polarized and unpolarized standard star measuremen ts were taken each ni ght in addition 
to the asteroid observations. Standards were drawn from lFossati et al.l (120071) as well as the 



standard hst published for Keck/LRISp0 which includes the Hubble standards ( Schmidt et al. 



19921 ). Standard measurements showed that induced instrumental polarization, systematic 
depolarization and instrumental crosstalk were all below the 0.1% level, and so are not 



^http: / /www2. keck. hawaii.edu/inst/lris/polarimeter/polarimeter. html 
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included in the data tables. See §[3] for further discussion on systematic error determination. 



m 



Our linear polariz at ion data are shown i n Figured] along with the Barb ara data presented 



Ce 



lino et al. 



(I2OO6I ). ICellino et al.l (120071 ) (these data also publi shed in iGil-Hutton et al. 



kooi ) as the Torino observations), and the CASPROF data from Icil-Hutton et al.l J2008h. 



In add ition, two comparison phase curves are also plotted. The first is the fit by lMuinonen et al 
( I2OO2I ) to polarization data from (24) Themis, a member of the B taxonomic class. B-type 
and C-type aste roids typically show the deepest negative polarization branches. The other 
curve is a fit by iGil-Hutton et al.l (120081 ) of data for (12) Victoria, a typical L-class object. 
Although different fitting models were used in these two cases, both models are qualitatively 
the same at the phase angles of interest. Not e that Aquitania was classified as an L-class 
and Barbara as an Ld-cla ss bv iBus Sz Binzell ((2002) whereas both are now L-class under 
the classification system of iDemed (120071 ) but show very different polarization features from 
other L-class objects. 



Instrumental Polarization Calibration 



Following the calibrations reported in iMasiero et al.l (120071 ) a quarterwave retarder was 
added in series with the halfwave retarder affecting the instrumental systematics, especially 
crosstalk between linear and circular polarization. To quantify these instrumental errors, ex- 
tensive lab benc h crosstalk calibratio ns were performed. Details of the calibration procedure 
can be found in lMasiero et al.l (120081 ). and are briefly summarized here. 



Figure [2] shows the measured polarization state for an input of pure linear polarization 
stepped through 360 degrees both without and with a fixed quarterwave plate in the light 
path (Figs [2^ and [2}3, respectively). Ideally the former test should show pure Q and U 
polarization as offset sinusoids peaking at 100% with zero V polarization, while the latter 
test should show offset Q and V (in this case, the U that is in the same phase as Q indi- 
cates a slight misalignment of the quarterwave plate at the input, while the out of phase 
U indicates crosstalk). A chi-squared minimizer was used to fit the measured variations in 
the Stokes vectors to determine crosstalk and depolarization. In all cases, the crosstalk and 
depolarization were found to be a few percent of the input polarization, so that for sources 
with "typical" polarizations (i.e. 5 — 10%) the errors are comparable to the desired statistical 
errors of ~ 0.1% polarization, as has been measured for polarized and unpolarized standards. 
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Discussion 



The ne w measuremen t s of (2 34) Barbara presented here agree with those from lCeUino et al, 

( 120061 ) and ICeUino et al.l (120071 ). Despite very different optical designs, data acquisition 



methods, observing circumstances, and physical locations, the consistency between our re- 



sults and those from previous work (e.g. ICellino et al.l (120061 ) ) indicate that comparisons and 



combinations of results from these instruments are legitimate. 

The similarity between the Barbara and (387) Aquitania polarization values clearly 
point toward Aquitania being a member of the new class of large inversion angle (LIA, or 
"Barbarian") as teroids that includes B arbara, (172) Baucis, (236) Honoria, (679) Pax and 
(980) Anacostia (iGil-Hutton et al.ll2008l). Barbara, Aquitan i a, and Anacostia all show strong 
2/i spinel features in their IR spectra ( Sunshine et al. 2007 : Burbine et al.lE992l ) implying a 
mineralogical origin to the unusual polarization properties of these asteroids. The other three 
LIA objects do not currently have publishe d IR spectral coverage at 2/x, but we predict that 
they too will show the same spinel feature. ISunshine et al.l (120071 ) explain that a spinel-rich 
spectrum can be used as a tracer of Calcium- Aluminum- Rich Inclusions (CAIs), one of the 
oldes t known materi a ls in t he Solar System, as determined from meteorite chemical analysis 



[e.g. 



Sunshine et all (120081 1). 



Gil-Hutton et al.l (120081 ) propose that the surfaces of these bodies may be coarse regoliths 



of a dark matrix mixed with smaller white inclusions. This mixing of two components with 
different albedos can alter the behavior of bulk polarization as a function of phase and explain 
the large inversion angle seen in these objects. The authors point out that their sample of 
LIA asteroids covers a wide range of semimajor axis-space (2.38 — 2.80 AU) and thus cannot 
be fragments from a single parent body. 

Table [2] s hows that Aquitania h as orbital and physical properties fairly s imilar to those 



of An acostia (iGil-Hutton et al.ll2008l ). This similarity was also pointed out by lBurbine et al. 



(119921 ) when they were both identified as spinel-bearing asteroids. However, it is very unlikely 
that these two objects originate from the same parent body, based on the velocity spread 
between their o r bits compared to tho se seen for other disrupted bodies and dynamical families 
( IWillmanll2008l : IZappala et al.lll995l ). In other words, if both Aquitania and Anacostia were 
fragments from the coUisional disruption of a common parent body, the difference in orbital 
elements would imply unrealistic ejection velocities of the order of many km/sec. Moreover, 
given the non- negligible sizes of Aquitania and Anacostia, such a collisional event should 
also be expected to have generated an important dynamical family that should be possible 
to identify even after a very long time. Such a family, however, is simply not found. 



It is worth noting that of the six asteroids now identified as LIA objects four have periods 
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longer than 20 hours when the average for these sizes is ~ 10 — 15 hours. Though still a 
victim of small number statistics, if this trend of long rotation periods holds for other LIA 
asteroids it may indicate that they are characterized by rotation p eriods significantly longer 
than the typical values found for most large Main Belt asteroids ( iPravec et al.ll2002l ). This 
would imply a highly porous, "fluffy" internal structure that is very efficient at absorbing 
energy from impacts without transferring it into an increase in rotati on rate. This kind o f 
structure is similar to what was observed for the Tagish Lake meteorite (jZolensky et al.ll2002l ) 
and what is expected for carbo naceous chondrites in which all CAIs found th us far have been 



observed (iBurbine et al.ll2002l ). This would support the theory proposed by ISunshine et al. 



(120081 ) that spinel-bearing asteroids, and thus probably all LIAs, are composed of pristine 
material from the beginning of the solar system. We note also that long spin periods might 
also be diagnostic of binarity, although any physical reason for a correlation between binarity 
and anomalous polarimetric properties remains unknown. 
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Table 1. Asteroid Polarization Measurements 



Target 


UT Obs Date 


V mag 


Tea:p (S6c) 


a 


Linear % Por 




Giro % PoP 


387 Aquitania 


2008-01-17 


12.7 


360 


21.2° 


-0.97 ±0.06 


91.36 ±2.57 


0.03 ±0.02 




2008-03-12 


11.7 


72 


15.6° 


-1.42 ±0.06 


88.95 ± 1.84 


-0.07 ±0.04 




2008-05-14 


11.3 


72 


17.1° 


-1.31 ±0.06 


89.58 ± 1.35 


0.00 ±0.03 




2008-06-11 


11.7 


60 


23.3° 


-0.82 ±0.07 


90.16 ±2.54 


-0.03 ±0.04 


234 Barbara 


2008-01-17 


14.3 


1620 


19.9° 


-1.24 ±0.06 


92.11 ±2.80 


0.06 ±0.03 




2008-03-12 


13.6 


720 


23.6° 


-0.93 ±0.04 


89.33 ± 1.61 


0.08 ±0.03 




2008-05-14 


12.2 


270 


13.7° 


-1.53 ±0.09 


87.98 ±2.70 


-0.01 ±0.03 




2008-06-11 


12.0 


108 


13.3° 


-1.67±0.13 


89.07 ±2.45 


0.06 ±0.03 




quoted errors are Icr 


statistical errors; 


systematic 


errors are 0.05%. 







Table 2. Orbital and Physical Parameters 



Asteroid 


a (AU)"^ 




i (deg)° 


H (mag)^ 


Diameter (km)^ 


Rot. Period (h)^ 


387 Aquitania 


2.739 


0.237 


18.14 


7.41 


100.51 


24.144 


980 Anacostia 


2.743 


0.200 


15.90 


7.85 


86.19 


20.117 



°- from MPCORB: hUp://www.cfa.h arvard.edu/iau/ MPCORB. html . ^ fromlTedesco et al.l hooj) . " from 
I http:/ / cfa-ww w. harvard. edu /i(m/lists/LightcurveDat.html\ 
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Fig. 1. — Polarization measurements of (387) Aquitania and (234) Barbara from this 
work (crosses and filled triangles, respectiv ely) compared to th e V and R ba n d Ba rbara 
measurements (open an d filled circles) from ICellino et al.l (120061 ) , ICellino et al.l (120071 ) and 



Gil-Hutton et al.l (l2008l ). The dashed and dotted lines show typical polarization-phase 
curves for B-type (e.g. (24) Themis) and L-ty pe (e.g. (12) Victoria) asteroids, respectively 



(IMuinonen et al.ll2002l : iGil- Button et all 120081 ). 
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Fig. 2. — Polarization crosstalk measurements for DBIP for two cases: (a) for input of 
pure linear polarization rotated through 360 degrees, and (b) for input of rotated linear 
polarization passed through a fixed quarterwave plate. Sinusoids for both cases were fitted 
using a chi-squared m inimizer, and are labeled with the Stokes vectors they represent. From 
Masiero et all J2008h 



